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Nonequilibrium transport via spin-induced sub-gap states in 
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We study low-temperature transport through a Coulomb blockaded quantum dot (QD) contacted 
by a normal (N), and a superconducting (S) electrode. Within an effective cotunneling model the 
conduction electron self energy is calculated to leading order in the cotunneling amplitudes and 
subsequently resummed to obtain the nonequilibrium T-matrix, from which we obtain the nonlinear 
cotunneling conductance. For even occupied dots the system can be conceived as an effective S/N- 
cotunnel junction with subgap transport mediated by Andreev reflections. The net spin of an odd 
occupied dot, however, leads to the formation of sub-gap resonances inside the superconducting gap 
which gives rise to a characteristic peak-dip structure in the differential conductance, as observed 
in recent experiments. 

PACS numbers: 73.63.Kv, 73.23.Hk, 74.45. -fc, 74.55. -|-v 



I. INTRODUCTION 

Magnetic impurities in normal metals are known 
to give rise to so-called Abrikosov-Suhl resonances, ^"21 
which in turn lead to the celebrated Kondo conductance 
anomaly observed in normal (N/barrier/N) tunnel junc- 
tions with magnetic impurities in the barrier,**'"^ as well 
as in normal (N/QD/N) cotunnel junctions based on 
Coulomb blockaded quantum dots (QD) holding an odd 
number of electrons .l^"^ Magnetic impurities in super- 
conducting metals, on the other hand, give rise to local- 
ized Yu-Shiba-Rusinov bound states inside the supercon- 
ducting gap,l^^lli2l which can be observed by measuring 
the local density of states (DOS) in scanning tunneling 
microscopy (STM) as sub-gap conductance peaks offset 
from the gap edge roughly by the magnitude of the ex- 
change coupling (cf. Refs. [3DH2S] and references therein) . 

In this paper, we explore the effects of spin-induced 
bound states in (S/QD/N) cotunnel junctions based on 
Coulomb blockaded quantum dots contacted to one su- 
perconducting and to one normal metal lead (cf . Fig. Ilk) . 
As we shall demonstrate, a coupling to the normal lead 
will broaden the localized bound states and the resulting 
scattering resonances will be reflected as characteristic 
sub-gap peaks, accompanied by pronounced dips at the 
gap edges in the nonlinear conductance. Basically, for 
a spinful quan tum dot, sub-gap transport via Andreev 
reflections^fiEZl probes the profile of the sub-gap states 
rather than simply the bare BCS DOS (cf. Fig. [ijj). 

Transport measurements on such S/QD/N systems 
in the cotunneling regime have bee n alr eady carried 
out.f^Stsni Most recently. Deacon et al^^ have indeed 
observed sub-gap conductance peaks for odd occupied 
S/InAs-QD/N devices, which they interpr et as signa- 
tures of Andreev energy levels inside the gap !^^ * '^^ ' Below, 
we argue that these peaks can be ascribed to Yu-Shiba- 
Rusinov resonances forming in a spinful cotunnel junc- 




FIG. 1: a) Sketch of the S/QD/N cotunnel junction com- 
prised by a (gray) quantum dot holding a net spin, which is 
tunnel coupled to an (orange) superconducting and a (green) 
normal metal lead, b) Illustration of the basic Andreev re- 
flection process giving rise to sub-gap transport and how it 
is enhanced by the presence of localized resonances induced 
by the spin on the quantum dot. The right (green) region 
illustrates the constant density of states (DOS) in the normal 
lead, shifted by the voltage V, and the left (orange) region 
illustrates the BCS DOS in the superconducting lead. The 
line is the local dot-electron DOS, with sub-gap resonances 
inside the gap and a pronounced dip at each gap edge. 



tion. This is consistent with the results of Refs. 1311321 
but allows for a simpler interpretation and calculation 
in terms of the Kondo model rather than the Anderson 
model. As we show, the experimental observation of en- 
hanced Andreev current in spinful dots, conductance dips 



near the gap edges and gate-dependence of the sub-gap 
peak positions can all be explained in terms of such spin- 
induced resonances. 

A number of works have addressed the problem of an 
Anderson impurity coupled to a single superconductor, 
either by numerical renormalization gr oup ( NRG) calcu- 
lationP^ESlor auxiliary boson methods j^^EZl and have ex- 
plored the intricate competition between Cooper-pairing 
and local correlations as a function of tunnel coupling F, 
charging energy U, and superconducting gap A. Even a 
simple no n-interacti ng (U — 0) model gives rise to sub- 
gap stateSjS^'^iMIMI which may affect the nonlinear con- 
ductance, and sub-gap states are thus sustained in many 
different parameter regimes and possibly with many dif- 
ferent characteristics. The present paper, however, fo- 
cuses on quantum dots in the cotunneling regime safely 
inside a Coulomb diamond where charge fluctuations are 
strongly suppressed. Restricting ourselves to the cotun- 
neling (Kondo) model, which is far simpler than the An- 
derson model, we retain crucial correlation features and, 
at the same time, we are able to capture the important 
physics of spin-induced bound states in a quantum dot 
setting, even out of equilibrium. 

As already mentioned, a finite coupling to the nor- 
mal metal lead will change the sub-gap bound states into 
broadened resonances but for low enough temperatures, 
Kondo correlations will become important and the res- 
onances will be either suppressed, or supplemented by 
a Kondo resonance pinned to the normal metal Fermi 
surface. The full S/QD/N problem is an inherently com- 
plicated problem,l3lBHB8] -^^iiiQij -^^g shall not attempt to 
solve here. In order to isolate and explore the observable 
consequences of the spin-induced sub-gap resonances, we 
neglect the log-singular terms arising from Kondo corre- 
lations with the normal lead, thus tacitly assuming the 
coupling between dot and normal lead to be sufficiently 
weak such that the corresponding Kondo temperature, 
Tk, is much smaller than than either temperature, T, or 
applied bias-voltage, V. Staying with an effective cotun- 
neling model, it would indeed be interesting to investi- 
gate the competition between these sub-gap resonances 
and Kondo instabilities in the regime A <|C Tk where 
nonlinear conductance has been reported^^ to be very 
different from the regime A ^ Tx, which we study here. 

After introducing the model, we explain how to get 
from lowest order self-energies to the current via the 
nonequilibrium T-matrix. Using this setup, we then 
start by investigating the case of a spinless even-occupied 
quantum dot, for which the effective cotunneling model 
takes the form of a simple potential scattering term. The 
nonlinear conductance is shown to be the same as for an 
ordinary S/N junction,'22MMl -^^hich is to say that the 
spinless quantum dot can be viewed as an effective S/N- 
cotunnel junction. 

Next, we treat the case studied originally by YrpSl, 
ShibcP^, and Rusinov^ of a classical spin, modelled by a 
spin- dependent potential scattering term and determine 
again the nonlinear conductance. In this case we can 
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FIG. 2: Comparison of the nonlinear conductance through 
an odd occupied quantum dot holding a net spin for different 
choices of coupling asymmetry. The spin exchange interaction 
is parametrized by J^^i = Jorar^'- A spin-induced bound 
state is seen as a sub-gap resonance for strong coupling to the 
superconducting lead rs S> tm (factor of 10 different) and for 
symmetric coupling, though broadened in that case. A dip 
appears at the superconducting gap, where one would expect 
a square root divergence like in the case of rjv ^ rg, where 
the normal lead is stronger coupled. 



still obtain exact analytical expressions for the current 
and comparing the expressions with the potential scat- 
tering case we find that already the classical spin leads 
to a peak below the gap accompanied by a square root 
dip at the gap edge (as opposed to the usual BCS square 
root divergence). 

In contrast to the case of potential scattering and the 
classical spin approximation, which can be solved exactly, 
the full quantum mechanical spin is investigated by T- 
matrix resummed perturbation theory within leading or- 
der in the cotunneling amplitude. The numerically de- 
termined results for the nonlinear conductance are sum- 
marized in Fig. [2] for three different regimes of coupling 
asymmetry. For stronger coupling to the S-lead, reso- 
nances similar to the Yu-Shiba-Rusinov states below the 
gap are very sharp and the dip at the gap edge well de- 
fined. As the coupling to the normal lead increases, these 
resonances become broader and start fiUing-in the dip. 
Eventually the sub-gap resonances merge with the dip 
and reproduce the usual BCS-profile at the gap edge. 



II. MODEL 

A. Model system 

We consider a quantum dot coupled to a normal and a 
superconducting lead (cf. Fig. Ilk). The highest partially 
occupied orbital on the dot is represented by a single- 



orbital Anderson model: 



H = Hs 



H 



N 



Ht + H, 



Dot, 



(1) 



with 



+ X!(^"C"kTCa-H + ^acJ^-HcLt) ' 

k 

akcr 

Hoot = ^ £ddid^ + Uud^ridi, (2) 

(7 

where a = N, S labels respectively a normal metal elec- 
trode (A^r — 0) and a superconducting lead with an or- 
dinary s-wave BCS DOS and a gap which is assumed to 
be real (Ag = A). This can be safely assumed, since 
the phase of the superconductor does not play a role 
in an S/N junction. Dot electrons of spin a are cre- 
ated by dj. in an orbital of energy Sd and with a mutual 
Coulomb interaction strength U. The tunnelling ampli- 
tude between dot and electrodes is denoted by t^, and 
the Coulomb blockaded dot is tuned by the gate-voltage 
{Vg ~ Ed) to hold a well-defined number of electrons. For 
a partial filling of this highest lying dot-orbital of one, 
i.e. a single electron on the dot, we thus assume that 
Ta = TTi^pltal"^ < max(-ed, Ed -I- U). 

In order to represent the exchange interaction between 
conduction electrons and the spinful quantum dot, it is 
necessary to augment the standard BCS Nambu-spinors 
to liberate the spin, from the charge. To this end, we 
introduce the four-spinors 



V'lk = ( 



t 



ak 



-akf *^a-kti <^ak^' Ca-k^ 



/ Cakt 

C^ 



(3) 

(4) 



k^ 



/ 



satisfying the anti-commutation relations 



{'0ak,7,l/'l'k'r,'} = '^"""^l 



kk'm 



,0 



/,(t) ,/,(t) 



{€Lr:^w} = -^-"^-kk'm 



vv 



(5) 
(6) 



where we have introduced the following set of 4 x 4 ma- 
trices: 



m° = 



10 
10 
10 
1, 

/ 
1 
10 
-10 



1 







10 
0-100 
10 
0-1 

/O 1 
10 
1 

Vo 1 0. 



(7) 



Within this notation, the Hamiltonian for the leads reads 
^" = jE^Uk"i"^k+j5]V'lA„m^Vk, (8) 



where a — S,N, using Ajv ~ and A5 = A. 

Since we focus entirely on cotunncling, we project out 
charge-flu ctuat ions by means of a Schrieffer- Wolff trans- 
formation.ESEBThis simplifies the Anderson model to the 
effective cotunneling model: 



^cotun — ^exch i -flo 



pot I 



(9) 



with 



H, 



(10) 



a']i.'ri^ .aWrji i—x^y^z 



where the 4-spinor notation has been supplemented by 
the following augmentation of the Pauli matrices: 



10 
0-1 
10 
0-100 



-i 
-i 
i 
i 



10 
0-100 
0-10 
1 



(11) 



The spin exchange interaction, Eq. ( 10 ), applies only to a 
spinful quantum dot holding a net spin, S = -^dlifrj'ada, 
whereas the potential scattering term 



i^, 



/,t 



pot = 2 E W^^'aV^a'k'^' 



1^'j^iPa]<i7] 



(12) 



applies to spinless quantum dots, as well as to spinful 
dots away from the particle-hole symmetric point in the 
middle of the relevant Coulomb blockade diamond. For 
a dot holding one electron, the two different cotunneling 
amplitudes are given by 



*J nn' — 



{ed + U){~edy 



and 



W„„, = 



(2ed + U)tatl, 
2{ed + U){~ed) 



(13) 



(14) 



where indeed Waa' = at the particle-hole symmetric 
point £d = —U/2. This result is readily generalized to 
any odd number of electrons on the dot, as long as only 
single-electron charge fluctuations are being eliminated 
by the Schrieffer- Wolff transformation. In the case of 
an even occupied dot, we denote the corresponding 2nd 
order cotunneling amplitude by W^^,. 



As an illustrative intermediate step, we shall also dis- 
cuss the case of a classical spin, where we replace the full 
exchange cotunneling term by a spin-dependent potential 
scattering term: 



and 



Tjcl. 

exch 



1 



Yl W^a'aV'l'k,,'K''7V'"k,, 



(15) 



Q'k'r;',ak77 



We stress that this is merely a simplification of H^xch 
corresponding to the limit of J — >■ and 5 — > oo, keeping 
the product, JS*, constant. This was in fact the problem 
considered e.g by ShibgP^ and as we shall demonstrate it 
already captures the essential physics of the spin-induced 
sub-gap resonances even for the full quantum spin. 



B. Unperturbed Green functions 

Since we want to study transport beyond the linear 
regime, we do the perturbation theory using Keldysh for- 
malism. We define the contour-ordered conduction elec- 
tron Green functions: 

Gakij.a'k'rj'lT,/) = -l (Tc | V'ak,, (T)V'l'k',,' ('^') } 

(16) 

In all expressions encountered below, the Green functions 
are summed over momentum. Therefore we start out 
by stating the unperturbed momentum summed Green 
functions for the leads. We have for the retarded and 
advanced Green functions 



k 



'TrvF0iD- \lu\ 



{uj±iO+)m° + Am'' 
' ^A2 - (w ± i0+)2 



(17) 



where D denotes the conduction electron half-bandwidth 
and we assume a constant DOS (pr. spin), vp. For the 
spectral function, we then have 



A(°)(a;) = 2tti'fO{D - \uj\)0{\uj\ - |A|)sign(a;) 

wm° -|- Am'' 
X — , 

y/uj^ - A2 



(18) 



Notice that the anomalous part of the DOS is an odd 
function of w, which falls off as A/w for large cu. For the 
normal lead, where A = 0, the spectral function reduces 
to 



A^°\uj) = 2t:vfO{D-\u\ 



(19) 



From the momentum-integrated spectral function, the 
lesser and greater Green functions are readily obtained 
as 



G(")^<(t.)=a(a;)/(a;), 



(20) 



G(")^>(c.)--a(c.)(l-/(c.)), 



(21) 



where / denotes the Fermi function. The voltage is ap- 
plied to the normal lead only, thus avoiding the compli- 
cation with a running phase in the superconducting lead, 
and therefore the voltage V enters only in the normal 
lead Green functions 

gK' = */(^ - m;^T/)m°„,yl(,°H^), (22) 

gK' - -*(1 - /(^ - m:;,T/))m°^,^(^")('^), (23) 

where the chemical potential shift has opposite sign for 
particle and hole components. 

C. Current and T-matrix 

Since we will be particularly interested in bound states 
(or resonances) as poles in the conduction electron T- 
matrix, it is convenient to express the current in terms 
of the T-matrix. To do this, we momentarily revert to 
the underlying Anderson model for which the current op- 
erator is found as the rate of change of the number of 
particles in the superconducting lead: 



dt{Qi 



(-e) 



2ih '^ 

k 



iJT>5k^"''^5k] 



(24) 



Introducing four-spinors for the dot electrons: 



<^t = ( 4, 



1^, a^. 



/'l\ 



V4/ 



(25) 



the tunneling term takes the following form: 



Q,k,T7 



(26) 



where the tunneling amplitude has been chosen to be 
real, which is always possible for a single-level model with 
one normal lead, since a phase can be absorbed by a 
gauge transformation. 

The expectation value of the current operator now 
involves the mixed 4x4 Nambu Green functions, 
i(T(7^5.^'(r)0jj(r')), and one can show that 



X ([G,,,,(^)G^?,,(^)]<- [4;U(^)G.;,,(^)]' 

(27) 

where Gd-,i~f is the dot electron Green function in spinor 
space. In Eq. (27) we use the shorthand [A{uj) B (uj)]'^ = 



A^{uj)B'^{ijj) + A'^{uj)B'^{uj) implied by Langreth rules 
and subsequent Fourier-transformation. From equations 
of motion, the dot electron Green function can be ob- 
tained from the conduction electron T-matrix as 



''a"^d;j;jj' l'"; '' ) ^'^r;r;'^^))'r;' -^ a;))?;' I'': '' J; 



(28) 



which may be inserted to obtain the following formula 
for the current: 



oo 



where 



^AfS,cff 



^ATMl ^ AT ^^ 






(34) 



(35) 



and 



T, 



SN 



{m^ + TisGf^^)i:^SN{m^ 






(36) 



where we have used the relation Q^ 



(29) 



.,,, = -(^<,,)^ The 
current now relies solely on the conduction electron T- 
matrix, which we calculate within the effective cotunnel- 
ing model derived from the Schrieffer- Wolff transforma- 
tion. Thus the details of the quantum dot enter only 
via the interaction with the leads and in the following 
we study Eq. ( 29 1 for Hpot and/or H^xch as defined in 
section III Al 

The T-matrix effectively sums up an infinite repetition 
of the one-particle irreducible self energy S, which we 
obtain either exactly (for spinless dots and for the limit 
of a classical spin) or to leading order perturbation theory 
in the exchange cotunneling amplitude. In general, the 
retarded or advanced SS component of the conduction 
electron T-matrix is found as 



nR/A, 



nr{uj) 



^R/A 
^SS,cff 



(^)(r 



G\ 



\'^)^%%{'^))' 



(30) 



with an effective 55* self energy which incorporates as 
well all processes going via the normal lead back into the 
superconducting lead: 



^R/A 



^R/A, 



"SS,cff(^) ~ ^SS i^) 



',{0),R/A 



+ I]^];-(c.)[[G^^'^'/-(.)]- - !]-/,-(..)] -^i]-/f(..) 



B./A, 



]-l^B./A, 



(31) 



The irreducible self energy Sqq,' depends on the details of 
the interacting system, thus determining the expressions 
for the T-matrix and the conductance as will be discussed 
now for the various cases. 



III. RESULTS - EXACTLY SOLVABLE 

A. Spinless dot - potential scattering 

For a spinless dot, as for example one of even occupa- 
tion with a singlet ground state, the effective cotunneling 
model has only the potential scattering term (12 1, and 



the irreducible conduction electron self energy is then 
exact already to first order in the potential scattering 
amplitude, W^^,: 






(37) 



Since the self energy is independent of frequency, it has 
S^ = E-^ = S^ and E< = 0, and the effective 55* Nambu 
matrix self energy therefore takes the following simple 
form: 



-,e,R/A 
■"SS,cff 



1 

■KVp 



(o-^ m° =F ial m°). 



(38) 



where =F refers to retarded and advanced components and 
the real dimensionless coefficients are defined by 



This effective self energy also enters the Dyson equation 
for the interacting SS Green function 



G 



RjA 
SS 



M = {[4°^'^/^H]-^-sf/^e«M}"- (32) 



with the definition 



'^Sn'^ns 



^ss 



'NN "t 1 



(39) 
(40) 



In the current formula Eq. ( 29 ) , the lesser T-matrix 



enters through the following combination of Nambu ma- 
trices (omitting the w-dependence in the following): 



1 QcCr 



-SS"- 



iO).A 
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TisG 



(o),< 
s 



rpR ^ 

^SS^ 



(0),< 
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^SS.cS 



^ss 



rrR 



T, 



SN 






SS,cff ■ 



^(0),-Ry.< 



G^ 



'N 



NS,eS 



^N ^NS,eS r^SS^ 



SS 
A 



(33) 



^^•^fW^^, 



(41) 



Here at is responsible for transmission and cr^ can be 
related to reflections back to the superconducting lead. 
Note that all expressions are in terms of the T-matrix 
Tss only. The T-matrix Tjvat cannot be written in the 
simple fashion of Eq. ( 38 ) since superconducting corre- 



lations are induced in the normal lead, whereas in T55 
the coupling to normal lead only modulates the already 



present anomalous contributions in the superconducting 
lead. 

If the normal lead was decoupled from the quantum 
dot, i.e. W^g — Wgj^ — 0, we would have no transport, 
i.e. a^ — 0, and a'f — ass- In the opposite limit where 
there are no reflections, i.e. Wgg — W^pf = 0, the sys- 
tem becomes equivalent to a tunnel junction, with u^ = 
and erf = ctns'^SN- In general, however, the cotunnel 
junction studied here involves both numbers, cr^,. Even 
though the self energy is exact in this case, it should be 
kept in mind that the validity of our initial Schrieffer- 
Wolff transformation relies on the fact that the four di- 
mensionless numbers, aaa' , are all much smaller than 
one. 

Inserting the above self energy (381 into Eqs. ( 33p6 ) 
and using current formula ( 29 1 , we can now obtain a 



closed analytical expression for the nonlinear conduc- 
tance at zero temperature. For voltages outside the gap, 
y > A, we have 



dl 
dV 



2e^ 4crf 



h V 



(42) 
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with the denominator defined as: 



FIG. 3: Current I ja for the spinless dot (even occupation) 
with equal dimensionless cotunneling amplitudes to S and 
N given by respectively a" = TvuFWy2 = 1,0.75,0.5,0.25, 
corresponding to a cotunnel junction transmission of a = 
0.8,0.56,0.2,0.015. Dotted line indicates the superconduct- 
ing gap V = A. 



V 



V>A 



= 2a^ + (1 + {atr + Kf) V^^WJW- (43) 



For voltages inside the gap, T^ < A, we have instead 

4e2 4(a|) 



dl 

dV 



h V 



(44) 



y<A 



with denominator 



n<A = (1 + {^•iY + i<YY{i - {v/^Y) 

+ A{atY{V/^Y- 



(45) 



In the limit of A — > or F 



oo. 



the differential con- 



ductance reaches the value ^a, with a cotunnel junction 



transmission 



4af 



(1 



^e\2 



+ {KY 



(46) 



which differs from that of a tunnel junctiorP^^Sl merely 
by the presence of the reflection terms Wgg and W^j^ 
comprising a'^ . Note that a rewriting of the full nonlinear 
conductance, ( 42]|45 ), in terms of this a makes it identical 
to the expression found in Ref. I^Sl for a tunnel junction, 
and to the one by BTK^ , when expressing their barrier 
parameter Z in terms of the transmission, Z^ = a~^ — 1. 
The conductance at F = A is exactly 4e^//i, whereas 
the zero-bias conductance is given by 



dl_ 

dV 



_ 4e2 
v=o h 



2cr? 



1 + «)2 + (a-) 



e\2 



(47) 



which agrees with the general result for the linear con- 
ductance of an N/S interface^*: G = Ae'^/h a^/(2 - a)^, 
valid for any transmission a. 
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FIG. 4: Differential conductance G = dl /dV for the even dot 
with same parameters as in Fig. [3] Inset: the nonlinear con- 
ductance Gs/N normalized by the value of Gjv/at — 2e^ /h a. 



Sub-gap transport for T^ < A is allowed by Andreev 
scattering processes, which proliferate with increasing co- 
tunneling amplitudes W^^ . This is clearly seen in Fig. [3] 
where we plot the TV curves for different tunneling am- 
plitudes. The corresponding nonlinear conductances are 
shown in Fig. |4] and seen to simply reflect the BCS DOS 
for the smallest chosen tunneling amplitude. 



B. Classical spin - spin dependent potential 
scattering 

We commence with an extension of the spinless dot, 
described by a potential scattering term W^, to a dot 



holding a classical spin, which we describe by Eq. ( 15 ) in 



terms of a spin-dependent potential scattering term W^ . 
In this case, the exact irreducible self energy is given by 






(48) 



In line with the spinless case, the effective SS self energy 
can be written as 



^s,R/A 



1 

'KVp 



{a^ m^ =F io-j m°) 



(49) 



where the same definitions (39) and (40) apply to the 



dimensionless coefficients, Ct/r, when simply replacing 
^aa' by W^aa'- ^otc that the spin symmetry is broken 
by m^ in contrast to the potential scattering case. 

As before, we can find a closed expression for the non- 
linear conductance at zero temperature. In Figs. [5] and [6] 
we first show the resulting I-V curves and corresponding 
nonlinear conductance for the same coupling strengths 
as used for the spinless case in Figs. |3]|4] We limit the 
plots to the case of symmetric couplings and return to in- 
vestigate asymmetric couplings for the full quantum me- 
chanical spin in the next section. For weak coupling, the 
conductance in Fig.[6]is similar to the potential scattering 
case and merely reflects the BCS density of states. For in- 
creasing coupling strength a sub-gap peak appears sym- 
metrically around zero bias and the conductance peak 
at V^ = A changes to a dip, as discussed below. As the 
peak moves closer to zero energy, it also becomes broader 
for our choice of symmetrically coupled leads. Note that 
for very strong coupling the two Yu-Shiba-Rusinov states 
mix and one sees only a broad peak centered at zero volt- 
age, which we should emphasize has nothing to do with 
a Kondo resonance. All of these features can be identi- 
fied in the following analytical formulas for the nonlinear 
conductance. 

As in the case of potential scattering we find for the 
classical spin exact expressions for the differential con- 
ductance from perturbation theory. For voltages outside 
the gap, t/ > A, we have 
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given by (43) with a'^ 
-^QorV 



replaced by a^ . In 



->■ cx) we find again that the con- 
ductance is given by 2e^//i a with the cotunnel junction 
transmission, Eq. (46), in terms of ct^^.. However, due to 



the term with a^ in the denominator, the differential con- 
ductance here cannot be expressed in terms of a alone. 

For tr^ « 0, i.e. weak coupling to the superconducting 
lead (Jss ~ Oj the conductance outside the gap is iden- 
tical to that in the spinless case with a BCS-like square 




FIG. 5: Current I /ct for a symmetrically coupled quantum 
dot occupied by an odd number of electron where the spin is 
treated classically. The coupling strengths are chosen g" = 
■rufW''/2 = 1,0.75,0.5,0.25 similarly to Fig.|3] 



root singularity at \eV\ ~ A. This changes rapidly with 
increasing a^, which changes the singularity at \eV\ = A 
into a dip with a square root increase with voltage (cf. 
Figs. [5 and|6]). At F = A the differential conductance 
takes the value 



dl 
dV 



4e^ 
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y=A h (a,«)2 + (a|) 



rS\2' 



(51) 



which reaches the 4e^//i, attained in the spinless case, 
when a^ <^ af, i.e. in the case where the dot is coupled 
much stronger to the normal, than to the superconduct- 
ing lead. On the contrary, for strong coupling to the 
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FIG. 6: Differential conductance G — dl/dV for the odd oc- 
cupied dot with a classical spin where the coupling strengths 
are the same as in Fig. |5| Inset: Conductance Gs/n normal- 
ized by Gn/n = 2e^//i a. 



superconducting lead, cr* :^ CTj , the value of the conduc- 
tance is suppressed and the dip is clearly visible. 

The spectral weight which has been removed at the 
gap edge can instead be found inside the gap. For volt- 
ages inside the gap, V < A, the conductance takes the 
following form: 






T, 



where 

Us = ±A 



(v-ujsr + riiv + ojsr + ry 



(l + (ai:)2-(a,^)2)2 + 4(aj:)2(af) 



(52) 

(53) 
(54) 



The sub-gap conductance shows two peaks of approxi- 
mately Lorentzian form, centered at energies iujs and 
having a width of Ts- In the case of vanishing coupling 
to the normal lead, i.e. af = and af. = a-gg, the reso- 
nances sharpen to form real bound states, located at: 



s^2 



Ws 



±A 



1 - JK) 
l + ia^rY' 



(55) 



This limit reproduces the case of a classical spin embed- 
ded in a bulk superconductor .liSHill^ Por weak coupling 
this bound state is thus offset from the superconduct- 
ing gap by roughly the interaction strength. For strong 
coupling the dependence changes, but interactions of the 
order of the band width will on the other hand conflict 
with the confinement of electrons on the quantum dot. 
The classical spin case is artificial in the sense that the 
spin symmetry is broken although the spin exchange in- 
teraction does not break this symmetry. In the next sec- 
tion it is shown that spin-induced sub-gap bound states, 
qualitatively similar to those found in this section, are 
still present if the spin is treated quantum mechanically. 



IV. RESULTS - QUANTUM SPIN 

A. Exchange cotunneling 

Whereas the calculation for the classical spin is ex- 
act, we have to rely on leading order perturbation theory 
when it comes to calculating the conduction electron self 
energy in the case of exchange cotunneling with a quan- 
tum mechanical spin. For zero magnetic field, the leading 
term in the conduction electron self energy is of second 
order in the exchange interaction: 



Sa,,.v(i - o =^ E E(^^w^'(^')) 



(0)S7™7'V 



in u 



E^«I"Gi'',7y(*-*V2"t- (56) 



Before presenting the general results, let us first take 
a look at the limit of negligible coupling to the normal 
lead, i.e. Jjvs ~ 0. With the experience from last section, 
we expect to find bound states showing up as zeros in 
the denominator of the T-matrix. These are located at 
energies uts for which 

0^det{[Gf^^^^{us)r~^%''{us)}, (57) 



and inserting the self energy ( 56 ) , we can identify bound 
states at energies 






(58) 



in terms of the dimensionless exchange cotunneling am- 
plitude 



9ss = -iT^vpJss ■ 



(59) 



Expression ([58| is also valid for the general case of a spin 
exchange interaction (10 1 with a higher spin (e.g. spin-1 
in an even dot) by replacing 3/4 with S{S -f 1). 

The bound state energies in Eq. (581 match those 
found for the classical spin ([55| (in the limit of vanish- 



ing coupling to the normal lead) if one replaces Wgg by 
\/3/4 J55. However, this analogy to the classical (or po- 
larized) spin case does not hold anymore when the cou- 
pling to the normal lead is non-zero. We shall return to 
this comparison in Section |IVB[ 

As for the classical spin, a finite coupling to the normal 
lead broadens the bound states into resonances. Unlike 
for the classical spin, however, we cannot give a closed 
analytical expression for the T-matrix nor for the non- 
linear conductance. To obtain the retarded T-matrix, we 
insert (56 1 into ((30]) and (31 1. The imaginary part of 



the T-matrix, which is proportional to the local DOS 
on the dot (cf. Eq. ([28|), is plotted in Fig. [TJ We 
parametrize the exchange cotunneling amplitudes from 
Eq. ( |l3l ) by Jaa' = JafaTa' , whcrc Jo = 45o/7rJ^F and 
fa = ia/tQ is the ratio between the tunneling amplitude 
to lead a = S,N and t^ = maxjtjv, is}- 

The spin-induced sub-gap resonance is seen to stay at 
the same position as long as the coupling to the supercon- 
ducting lead is stronger than the coupling to the normal 
lead, i.e. 1 = rs > r^v in Fig. [7^. The peak though gets 
broader and lower with increasing coupling, r^, to the 
normal lead. In Fig. [TJd we show the behavior when the 
spin is stronger coupled to the normal lead by reducing 
the coupling to the superconducting lead 1 = tat > rg. 
Besides a strong suppression of the overall value, we also 
observe that the resonance moves out towards the gap 
edge and, for very weak coupling to the superconduct- 
ing leads, eventually gives back spectral weight to re- 
construct the usual square root divergence at the super- 
conducting gap A. These characteristics will be present 
again in the nonlinear conductance through the quantum 
dot as discussed later on in this section. 
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FIG. 7: Imaginary part of the T-matrix, - ImlTTUFTss], for a 
spinful quantum dot asymmetrically coupled to a supercon- 
ducting and normal lead with go = -kvfJq/^ = 0.3 leading to 
ids ~ 0.8735A for rjv = 0; a) rg = 1 is kept constant and 
the coupling to the normal lead vn is changed; b) rjv = 1 is 
kept constant and the coupling to the superconducting lead 
is reduced. 



1. Gate dependence of the spin-induced sub-gap resonance 
energy 

At the particle-hole symmetric point, Sd = —U/2, the 
spinful dot derived from the Anderson model has no po- 
tential scattering term, i.e. W = at this point (cf. 
Eq. ([T4|). Leaving this point by adjusting the gate- 
voltage (which is proportional to Ed), however, the po- 
tential scattering term has to be included, and it is clear 
from Eqs. ([l3 14 1 that both J and W will in fact in- 
crease in magnitude until eventually the Schrieffer- Wolff 
transformation breaks down as one comes too close to 
a charge-degeneracy point for the Coulomb blockaded 
quantum dot. 

To investigate the dependence of the resonance fre- 
quency, COS, on gate- voltage, we again consider the case 



FIG. 8; Gate dependence of the spin-induced sub-gap reso- 
nance a;s/A in terms of a; = \ed\/U, and with gss = 0.3; red 
(full) line including potential scattering and blue (dashed) 
line without. Black (dotted) line indicates the energy of the 
superconducting gap. 



where the dot is coupled only to the superconducting 
lead. As before, this is found as a root for the denomi- 



nator in the T-matrix (cf. Eq. (57)), but now we have to 
include also the first order term from potential scattering, 
i;(i) = Wqq'TO°/2, in the irreducible self energy. Doing 
this, one finds that the T-matrix pole will be located at 



UJs 



±A, 



'{i-v^9lsr 



^ss 



(1 + 3/451^)2 +w;^ 



(60) 



ss 



where wss = t^i'fWss/'^- For wss — and S = 1/2, 
we thus reproduce Eq. ( [58| for a spin coupled to a su- 
perconductor and for gss = we find ws = ±A, i.e. no 
spin-induced sub-gap state if there is no exchange cou- 
pling and only potential scattering as discussed in sec- 
tion HIEl 

In Fig. [8]we show how the position of the sub-gap state 
changes inside the Coulomb diamond as a function of 
X = \ed\/U. As the coupling strength defined in Eq. (13) 
increases towards the edges of the diamond (i.e. a; = 
and a; = 1), the bound state moves closer to zero. At the 
same time we of course expect the width of the peak to 
increase both due to the added influence of the potential 
scattering term but first and foremost due to the increase 
in cotunneling amplitudes as one moves away from the 
particle-hole symmetric point {x = 1/2). Notice that the 
potential scattering term, with the gate dependence as 
defined in Eq. ( 14 ) , makes practically no difference until 
perturbation theory breaks down anyway. Interestingly, 
the gate dependence shown in Fig. [8] is very similar to 
that reported in recent experiments on N/QD/S^S and 
S/QD/S junctions.'52l 

For the rest of the paper we solely deal with the 
particle-hole symmetric point where the potential scat- 
tering is zero. 
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FIG. 9: Comparison of the differential conductance through 
a spinful quantum dot, treated quantum mechanically, taking 
into account all terms (red data points) or the the weak cou- 
pling term only (black line) for symmetric coupling rs/r^ = 1 
and go = 0.3 . Furthermore, conductance treating the spin 
classically, dashed (green) line, with a coupling of W^^i = 

V3/4 v;;^. 



B. Transport via spin-induced sub-gap resonances 

To derive the nonlinear conductance in the case of a 
quantum mechanical spin, we need to include the full 
Keldysh matrix structure in Eq. ( 29 1 and ( 33 1 . Therefore 



we cannot provide an analytic expression for this case, 
but the calculations are straightforward and do not need 
extensive numerical effort. 

We can distinguish between a hierarchy of contribu- 
tions in Eq. (33). Summing over ttiJ? in Eq. (29), 

T^gGg cancels out while ^ss eff^ss contributes to 
the current already to second order in the spin exchange 
interaction Jaa'- This is the dominant contribution in 
the conductance as illustrated in Fig. [9] There "weak 
coupling" refers to calculating the nonlinear conductance 
including '^ss eff^ss only, while "full" refers to includ- 
ing furthermore all contributions in Eq. ( 33 ) which are 
of fourth order in lowest order in Jaa' 



As can be seen 
in Fig. [9J the weak coupling term actually overestimates 
the height of the spin-induced sub-gap resonance and the 
value at the edge of the superconducting gap slightly. 
However, for the coupling go — 0.3 (cf. Fig. [9]) studied in 
the following, this provides a very good approximation 
and therefore Figs. [T0p3| are calculated with the con- 
tribution from ^ss eff^s O'^ly- ^^^ stronger coupling 
the deviation is more severe. Furthermore, Kondo cor- 
relation effects have to be taken into account for strong 
coupling to the normal lead and this regime is not dis- 
cussed here. 

Fig. [9] also shows a comparison with the conductance 
in the case of a spin treated classically. For symmetri- 



cally coupled junctions we find that Waa' — ■\/3/4 \/ Jaa' 
provides a reasonable agreement although the energy of 
the spin-induced sub-gap state is slightly shifted and 
the value of the conductance is in general overestimated. 
Whereas a linear relation, Wgg = •\/3/4 Jgg, is a good 
approximation for a spin decoupled from the normal 
lead, a square root dependence {Waa' — ■\/3/4-\/ Jaa' ) fits 
the quantum mechanical case for symmetrically coupled 
junctions. In practise, we can always find a value for the 
classical spin case, which fits rather well the full spin-flip 
scattering case, but (as already clear from the two simple 
limiting cases) there is no obvious systematics involved 
and the value is strongly dependent on asymmetry and 
strength of the coupling to the leads. Nevertheless, we 
claim that the classical spin case provides good quali- 
tative insight into the problem of a spin coupled to a 
superconductor. 

Since the conductance is directly related to the lo- 
cal DOS, i.e. Im[G55], it is not surprising, that sub-gap 
states are observed in the transport through an N/QD/S 
cotunnel junction. As already illustrated in Fig. [TJ the 
local DOS of the superconductor is probed by Andreev 
scattering processes to the normal lead. The enhanced 
spectral density at the energy of the spin-induced sub-gap 
state leads to a sub-gap peak in the differential conduc- 
tance and a reduced DOS at the superconducting gap A 
as reflected in dl /dV . Comparing Figs. 10 and 11 with 
Fig.Mwe find significant agreement, since Im[T|^^]ls (be- 
sides a prefactor proportional to the interaction) given by 
the interacting local DOS. 

As shown in Fig. [TO] the spin-induced sub-gap state 
stays at roughly the same energy ws given by Eq. ( 58 ) 
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FIG. 10: Differential conductance G — dl/dV for increasing 
coupling to normal lead, rs = 1 and go = 0.3. The conduc- 
tance is renormalized by the "transport" coupling strength, 
JffS^ across the dot in order to compare the different curves. 
Inset: Illustration of a spin state coupled to the supercon- 
ducting lead, thus forming a bound state, which is probed by 
the normal lead. 
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(for a spinful quantum dot decoupled from the normal 
lead) as long as the superconducting lead is stronger cou- 
pled to the dot than the normal lead, rs > r^f. This can 
be understood schematically as illustrated in the inset of 
Fig. 10 If the spin is coupled strongly to the supercon- 



ducting lead, a bound state is created and can be probed 
in transport through the N/S setup with a weakly cou- 
pled normal lead. 

On the contrary, if the normal lead is stronger cou- 
pled to the impurity, r^v 3> J's, we only probe the BCS 
superconducting DOS (i.e. no spin-induced sub-gap reso- 
nances). This is sketched in the inset of Fig. 11 Starting 
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FIG. 11: Renormalized nonlinear conductance G/Jj^s for 
increasing coupling to the superconducting lead and rjv = 
1; Inset: Illustration of a spin state strongly coupled to the 
normal lead and thus only the superconducting DOS is probed 
in transport. 

from a nonlinear conductance displaying a clear sub-gap 
peak for rjq = rs, the resonance moves closer towards 
the energy of the superconducting gap A as the coupling 
to the normal lead is increased, rjv > rs in Fig. 11 If the 



coupling to the normal lead dominates, no sub-gap res- 
onance can be distinguished from the square-root singu- 
larity at A. As was also illustrated in Fig. [2] and is shown 
again in Figs. [lO] and 111] an experiment of the transport 
through an N/QD/S junction can therefore have very 
different signatures depending on the asymmetry of the 
coupling. 



It is expected from the expression oi ujs, (58), that the 



spin-induced sub-gap bound state moves into the gap for 
increasing coupling strength. This is illustrated in Fig. |12| 
choosing symmetric coupling to the superconducting and 
normal lead. For go = 0.1 the sub-gap resonance state 
is present, but hidden at V = A. However, already for 
go — 0.2 the conductance is seen to change into a peak 
at cjjs ~ 0.94A and an associated dip at A instead of 
the square root divergence of the clean superconducting 
DOS. Note that the conductance value at A decreases 
with increasing coupling to the superconducting lead. 
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FIG. 12: Renormalized differential conductance G/Jf^g for 
increasing coupling strength, go = 0.1,0.2,0.3,0.4, and 0.5, 
where both leads are symmetrically coupled, vn ~ rs ~ 1- 




FIG. 13: Renormalized current I / J%g for the same parame- 
ters as Fig. [12] 



Since we have chosen symmetric coupling in Fig. |12[ the 
increasing coupling to the normal lead causes a concomi- 
tant life time broadening of the peak inside the gap sim- 
ilar to the case of a classical spin. 

Finally, Fig. [13] shows the current for the same pa- 
rameters as the conductance in Fig. [T2j The dip in the 
conductance at the superconducting gap is hardly vis- 
ible in the current, where it should show up as a kink. 
Most interestingly, the current has a sharp increase at the 
energy of the spin- induced sub-gap resonance W5. This 
could easily be misinterpreted as a reduced superconduct- 
ing gap A, since there is no obvious difference between 
the line shapes of the curves with go =0.1 and go = 0.4 
in Fig. [13] Therefore, a clear distinction between a re- 
duced superconducting gap and a spin-induced sub-gap 
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resonance signature is best seen in the differential con- 
ductance. 



V. SUMMARY AND DISCUSSION 

We have investigated the transport characteristics of a 
Coulomb blockaded quantum dot sandwiched between a 
normal and a superconducting lead. The focus has been 
on the difference between a dot with an even number of 
electrons and one with an odd number of electrons or, 
in more general terms, the difference between dots with 
zero or finite spin. 

We have restricted the calculations to the cotunnel- 
ing regime, where charge fluctuations are strongly sup- 
pressed. The effective model is in this situation a "co- 
tunnel junction", where for even occupancy (or a spin- 
less dot) one obtains a simple tunneling Hamiltonian for 
tunneling between N and S, plus reflection terms (N to 
N, and S to S). In the language of the SchriefFer- Wolff 
transformation, this is known as potential scattering, see 
Eq. (12 1. In contrast, for odd occupancy (or a spinful 



dot) there is an additional term, namely the Kondo, or 
exchange-cotunneling, Hamiltonian, see Eq. (10), where 



the tunneling electrons couple to the spin on the dot and 
hence may induce spin flips. Furthermore, at the particle- 
hole symmetric point, i.e. in the middle of the odd occu- 
pancy diamonds, the potential scattering term is absent 
and only the Kondo Hamiltonian remains. 

Because the effective even-occupied dot Hamiltonian 
is quadratic, the current-voltage characteristic can be 
calculated exactly, and a general expression that only 
depends on the transmission (i.e. the normal state con- 
ductance) has been derived here and also pre viously in 
the literature for N/S tunnel-junctions.'^ZESIMI Por weak 
transparency, the resulting differential conductance is 
suppressed inside the gap and develops the well-known 
BCS square-root singular DOS at the gap edge. With 
larger transparency the conductance is finite inside the 
gap due to Andreev reflections, in full accordance with 
the BTK model.27 

The problem with a spinful odd-occupied dot Hamil- 
tonian cannot be solved analytically, and in order to 
study the influence of the quantum spin we calculate 
the lead electron self energy to second order in pertur- 
bation theory which is subsequently summed in the T- 
matrix. For this purpose, a general expression for the 
current in terms of the lead electron T-matrix has been 
derived. For weak tunnel couplings (also weak enough 
that Kondo physics is not relevant, TV ^ ^), this ap- 
proximation captures the important physics. Interest- 
ingly, we show that the calculation gives results which 
are qualitatively similar to a "classical" approximation, 
where the spin operator is replaced by a static magnetic 
moment. The static spin approximation could for exam- 
ple result from a mean-field approach with an (unjusti- 
fied) breaking of spin- rotational symmetry. Nevertheless, 
the classical spin model provides good insight by analogy 



to the well-known Yu-Shiba-Rusinov bound states, which 
appear when electrons in a superconductor scatter off a 
magnetic impurity. In the classical case, the bound state 
energy can be determined exactly, and is located inside 
the gap. For the quantum spin case there is also a bound 
state inside the gap, but it does not rely on the unphysical 
assumption of broken spin rotation symmetry. One may 
understand the origin of this spin- induced sub-gap bound 
state from dynamically reduced superconducting correla- 
tions due to the presence of an uncompensated magnetic 
spin. A Cooper pair consisting of a spin-up and a spin- 
down electron is roughly speaking both repelled and at- 
tracted by the impurity. Thus when breaking a Cooper 
pair, the energy of the localized excited state is smaller 
than the superconducting gap by roughly the exchange 
energy as a quasiparticle can gain energy by a spin-flip 
process. 

The spin-induced resonance states have profound con- 
sequences for the transport characteristics, because they 
give rise to a sub-gap feature in the differential conduc- 
tance. This feature moves further inside the gap for 
stronger coupling to the superconducting lead. The posi- 
tion of the resonance depends only weakly on the normal 
lead coupling, which on the other hand serves to broaden 
the resonance. In measurements, the sub-gap conduc- 
tance peak could easily be mistaken for a reduced-gap 
peak, and even more so since at eV = A there is a dip 
instead of a peak. This characteristic dip-peak structure 
has already been seen in experiments in Refs. I30|52l An- 
other clear prediction resulting from the calculation is 
that the position of the sub-gap structure moves away 
from the gap edge, to lower voltages, when the gate po- 
tential is tuned away from the particle-hole symmetric 
point in the middle of the diamond. 

The notion of spin-induced sub-gap resonances in 
the differential conductance of S/QD/N junctions leaves 
some interesting questions unanswered. First of all, it is 
not understood how the Yu-Shiba-Rusinov state gradu- 
ally changes into a Kondo resonance with increased cou- 
pling to the normal lead. Secondly, the influence of an 
applied magnetic fleld is not clear. Naively one would 
expect the sub-gap conductance peaks to split in a B- 
fleld. However, from the classical spin model we know 
that a spin dependent cotunneling model gives a very 
similar IV curve and from this analogy one expects lit- 
tle dependence on magnetic field, other than the overall 
suppression of superconductivity, of course. Finally, we 
mention the interesting problem of a spinful QD coupled 
to two superconducting leads, where the interplay be- 
tween spin-induced sub-gap resonances and multiple An- 
dreev refiections c an be expected o give rise to unusual 
transport features.'^SEll 
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